Geometric phases and criticality in spin chain systems 
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A relation between geometric phases and criticality of spin chains is established. As a result, we 
show how geometric phases can be exploited as a tool to detect regions of criticality without having to 
undergo a quantum phase transition. We analytically evaluate the geometric phase that correspond 
to the ground and excited states of the anisotropic XY model in the presence of a transverse 
magnetic field when the direction of the anisotropy is adiabatically rotated. It is demonstrated 
that the resulting phase is resilient against the main sources of errors. A physical realization with 
ultra-cold atoms in optical lattices is presented. 
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Since the discovery by Berry geometric phases in 
quantum mechanics have been the subject of a variety 
of theoretical and experimental investigations Possi- 
ble applications range from optics and molecular physics 
to fundamental quantum mechanics and quantum com- 
putation 0. In condensed matter physics a variety of 
phenomena have been understood as a manifestation of 
topological or geometric phases 0, HJ El 0- An in- 
teresting open question is whether the geometric phases 
can be used to investigate the physics and the behavior of 
condensed matter systems. Here we show how to exploit 
the geometric phase as an essential tool to reveal quan- 
tum critical phenomena in many-body quantum systems. 
Indeed, quantum phase transitions are accompanied by a 
qualitative change in the nature of classical correlations 
and their description is clearly one of the major inter- 
ests in condensed matter physics 0, 0. Such drastic 
changes in the properties of ground states are often due 
to the presence of points of degeneracy and are reflected 
in the geometry of the Hilbert space of the system. The 
geometric phase, which is a measure of the curvature of 
the Hilbert space, is able to capture them, thereby reveal- 
ing critical behavior. This provides the means to detect, 
not only theoretically, but also experimentally the pres- 
ence of criticality without having to undergo a quantum 
phase transition. 

In this letter we analyze the XY spin chain model and 
the geometric phase that corresponds to the XX criti- 
cality. Since the XY model is exactly solvable and still 
presents a rich structure it offers a benchmark to test 
the properties of geometric phases in the proximity of 
a quantum phase transition. Indeed, we observe that, 
an excitation of the model obtains a non-trivial Berry 
phase if and only if it circulates a region of criticality. 
The generation of this phase can be traced down to the 
presence of a conical intersection of the energy levels lo- 
cated at the XX criticality. This geometric interpretation 
reveals a relation between the critical exponents of the 
model. The insights provided here shed light into the 
understanding of more general systems, where analytic 



solutions might not be available. A physical implemen- 
tation is proposed with ultra-cold atoms superposed by 
optical lattices [ill It2T ] . It utilizes Raman activated tun- 
neling transitions as well as coherent drive of the system 
via Bragg scattering. The independence of the generated 
phase from the number of atoms, its topological nature 
and its resilience against control errors makes the pro- 
posal appealing for experimental realization. 

Consider the one dimensional spin- 1/2 XY model, with 
N spins, in a transverse magnetic field. This chain has 
nearest neighbor interactions with Hamiltonian given by 
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where M = (N — l)/2 for N odd. In particular, we are 
interested in the Hamiltonian that can be obtained by 
applying a rotation of <fi, around the z-direction, to each 
spin 



M 
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The family of Hamiltonians that is parameterized by <j> 
is clearly isospectral and, therefore, the critical behavior 
is independent from <f>. This is reflected in the symmet- 
ric structure of the regions of criticality shown in Fig- 
ure n I n addition, due to its bilinear form, H (<ft) is 
7r-periodic in <f>. The Hamiltonian H((f>) can be diag- 
onalized by a standard procedure, which can be sum- 
marized in the following three steps: (i) the Jordan- 
Wigner transformation, which converts the spin oper- 
ators into fermionic operators via the relations, a/ = 
(rim</ a m)( a i + *°f)/2; (ii) their Fourier transforma- 
tion, d k = ^Ez^e-* 2 -"^, with k = -M,...,M; 
and (iii) the Bogoliubov transformation, which defines 
the fermionic operators, bk — dk cos % — id_ k e 2u ^ sin 
where the angle Ok is defined by cosdk = efc/Afc with 
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FIG. 1: (a) The regions of criticality are presented as a func- 
tion of the Hamiltonian parameters A, 7 cos 20 and 7 sin 20. 
The geometric phase corresponding to the ground state (b) 
and the relative one between the ground and excited state (c) 
as a function of the path parameters A and 7. Values of the 
geometric phase corresponding to the loops Ci , C2 and Cs in 
(a) are also indicated. 



procedures diagonalize the Hamiltonian to a form 

M 

H(<t>)= A ^ k h- (2) 

k=-M 

The ground state of H(cf>) is the vacuum \g) of the 
fcrmionic modes, bk, given by 



i-9> = n ( c 
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where |0) fc and |l) fc are the vacuum and single excitation 
of the k-th mode, dk, respectively. The regions of crit- 
icality which appear when the ground and first excited 
states become degenerate are shown for this model in Fig- 
ure ^a). The XX model, which corresponds to 7 = 0, 
has criticality region along the line between A = 1 and 
A = — 1. The region of the Ising model phase transition 
are the two planes at A = 1 and A = — 1. The interesting 
paths of evolution for generating a Berry phase are those 
for which the state of the system can evolve around a 
region of criticality. Here we shall focus on the criticality 
region corresponding to the XX model 1 which can be 
encircled by adiabatically varying the angle <f> from to 
it. The corresponding Berry phase of ground and first 
excited states can be evaluated as a function of A and 7. 

Using the standard formula it is easy to show that the 
Berry phase of the ground state \g) is given by 



h)- (4) 



This result can be understood by considering the form of 
\g), which is a tensor product of states, each lying in the 



two dimensional Hilbert space spanned by |0) fc \0)_ k and 
For each value of k(> 0), the state in each of 
these two-dimensional Hilbert spaces can be represented 
as a Bloch vector with coordinates (20, 9k). A change 
in the parameter determines a rotation of each Bloch 
vector about the z direction. A closed circle will, there- 
fore, produce an overall phase given by the sum of the 
individual phases as given in (jJJ and illustrated in Fig- 
ure ^b). It is worth noticing that when 7 is small, the 
overall phase is strongly dependent on the exact values 
of A and 7 due to the (large) proportionality factor M. 
On the other hand, when 7 > A, each pair of k and — k 
states from (0 contributes almost equally to the over- 
all Berry phase (see Figure^b)). In this case, we obtain 
<p g /M « 7r which makes effectively the dependence of the 
overall Berry phase on the parameters 7, A much weaker. 

Of particular interest is the relative geometric phase 
between the first excited and ground states given by 
the difference of the Berry phases acquired by these two 
states. The first excited state is given by |efe ) = b\ \g), 
i.e. it is the single excitation of the b fermionic field, with 
ko corresponding to the minimum value of the energy A&. 
The relative geometric phase then becomes 



Peg = <Pe ~ fg = 



(5) 



If the parameter is adiabatically changed from to tt, 
we obtain ip eg = — 7r(l — cos#fc ) which, in the thermo- 
dynamical limit (N — ► 00), takes the form 
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for |A| > 1 -7 2 
for |A| < 1 - 7 2 
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As can be seen from Figure ^c), the most interesting 
behavior of ip eg is obtained in the case of 7 small com- 
pared to A. In this case ip eg behaves as a step function, 
giving either it or phase, depending on whether |A| < 1 
or |A| > 1, respectively. This behavior is precisely de- 
termined from whether the corresponding loop encloses 
a critical point or not and can be used as a witness of its 
presence. In particular, in the |A| < 1 — -f 2 case the first 
term corresponds to a purely topological phase, while the 
second is a geometric contribution [3(. 

This is related to a more general property: a non-trivial 
Berry phase is generated when the closed loop in the pa- 
rameter space is spanned in the vicinity of a degeneracy 
point. The topological nature of this phase is evident 
in the case of a loop contracting to a point. If, in this 
limit, the Berry phase remains non- vanishing, then the 
contraction point identifies the position of the degener- 
acy. As a critical point is also a point of degeneracy it is 
apparent that this procedure can be used to locate the 
criticality of an Hamiltonian. 

To better understand the properties of the relative ge- 
ometric phase in our model, let's focus on the region of 
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parameters with 7 <C 1, i.e. for a loop around the region 
of criticality which contracting to a point. In this case, 
it can be shown that the Hamiltonian, when restricted 
to its lowest energy modes dk and d-k , can be casted 
in a real form and, for |A| < 1, its eigenvalues present a 
conical intersection centered at 7 = 0. It is well known 
that when a closed path is taken within the real domain 
of a Hamiltonian, a topological phase shift tt occurs only 
when a conical intersection is enclosed. In the present 
case, the conical intersection corresponds to a point of 
degeneracy where the XX criticality occurs and it is re- 
vealed by the topological term in the relative geometric 
phase tp eg . It is worth noticing that the presence of a 
conical intersection indicates that the energy gap scales 
linearly with respect to the coupling 7 when approach- 
ing the degeneracy point. This implies that the critical 
exponents of the energy, z, and of the correlation length, 
v, satisfy the relation zv = 1 which is indeed the case for 
the XX criticality [Hj. 

From an experimental point of view, this method of 
obtaining a geometric phase is robust against relevant 
sources of errors. The most significant one originates 
from undesired transitions induced by the time depen- 
dence of the Hamiltonian. In any adiabatic process the 
occurrence of these transitions is hindered primarily by 
the presence of a finite energy separation between neigh- 
boring levels. However, as in many condensed matter 
systems, these energy separations may decrease quite 
rapidly as the thermodynamical limit is approached. In 
fact, in this model, the gaps between the energies of the 
lowest excited states scale typically as 7V~ 2 , which might 
substantially reduce the applicability of the adiabatic ap- 
proximation in a truly many-body setup. However, in the 
adiabatic evolution proposed here, all the transitions be- 
tween closely separated energy levels are prohibited due 
to symmetry reasons. It is apparent form equation 
that the adiabatic evolution is generated by the Hamil- 
tonian term H oc J2k a k ~ J2k dj.dk, which preserves the 
excitations in each fcrmionic mode, d^. In fact, the ex- 
citations I ek) = b\ \g) of model Q are also eigenstates 
of H . Therefore, the only allowed transitions induced by 
H are between the ground state and the doubly excited 
states frj^tfc |<?)- The corresponding finite energy gap 2A^ 
is then sufficient to adiabatically prevent these unwanted 
transitions, even in the thermodynamical limit. More- 
over, this symmetry constraint guarantees the viability 
of the adiabatic approximation also when energy level 
crossings occur between singly excited states |efc)'s. 

Another important feature of this model is the behav- 
ior of the phase under errors in the initial preparation. 
The smooth dependence of the geometric phase on the 
momentum x = ^irk/N, insures that a small deviation 
$ = X ~~ Xo from the value \o — /N (corresponding 
to first excited state, i.e. the minimum of the function 
A(x)) does not drastically affect the interference pattern. 
Indeed a perturbation analysis shows that states with a 



small momentum difference, <5, from xo acquire geometric 
phases with values shifted by an amount proportionally 
to 8. In particular, due to the independence of the func- 
tion A(%) from the system size, N, the spread of the ge- 
ometric phase only depends on 6 and remains unchanged 
in the thermodynamical limit. Naturally, this effect is 
accompanied by a reduction of the visibility and a shift 
in the value of the measured phase, which, however, only 
depend on the preparation and measurement accuracies 
and does not depend on the system size. 

This construction, apart from its theoretical interest, 
offers a possible experimental method to detect criti- 
cal regions without the need to cross them, which un- 
dermines the ability to identify the state of the sys- 
tem due to the presence of degeneracy. In particu- 
lar, we shall implement this model with optical lattices. 
To this end, consider two bosonic species labelled by 
a = a,b that can be given by two hyperfine levels of 
an atom. Each one can be trapped by employing two 
in-phase one dimensional optical lattices. The tunnel- 
ing of atoms between neighboring sites is described by 
V = -J2ia( J ° a L a (i+i)'? + H.c.). When two or more 
atoms are present in the same site, they experience colli- 
sions given by = %^aL a L' a ^' a ^- We sna11 
consider the limit J <^U where the system is in the Mott 
insulator regime [l^j ] with one atom per lattice site. In 
this regime, the effective evolution is obtained by adia- 
batic elimination of the states with a population of two or 
more atoms per site, which are energetically unfavorable. 
Hence, to describe the Hilbert space of interest, we can 
employ the pseudospin basis of ||) = \nf = l,n| = 0) and 
||) = \nf = 0, tk = 1), for lattice site I, and the effective 
evolution can be expressed in terms of the corresponding 
Pauli operators. It is easily verified that when the tun- 
neling coupling of both atomic species is activated, the 
following exchan ge i nteraction is realized between neigh- 
boring sites [niii nit. 

Hi = _^^(afaf +1+ afaf +1 ). (7) 

Uab { 

In order to create an anisotropy between the x and y spin 
directions, we activate a tunneling by means of Raman 
couplings [T3.ll5| . Application of two standing lasers L\ 
and L2, with zeros of their intensities at the lattice sites 
and with phase difference <j), can induce tunneling of the 
state |+> = (e-^/ 2 |a) + e^/ 2 |6))/\/2. The resulting tun- 
neling term is given by V c — J c J2i( c l c i+i + H.c), where 
ci is the annihilation operator of |+) state particles. The 
tunneling coupling, J c , is given by the potential barrier 
of the initial optical lattice superposed by the potential 
reduction due to the Raman transition. The resulting 
evolution is dominated by an effective Hamiltonian given, 
up to a readily compensated Zeeman term, by 

^ U ab 
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where g((f>) was defined in £[J. Combining the rota- 
tionally invariant Heisenberg interaction H\ with Hi 
gives the rotated XY Hamiltonian described by equa- 
tion ||TJ, where the parameter 7 is given by J^/(2eU a b) 
and e = (2J a Jb + J^/2)/U a b is the overall energy scale 
multiplying the Hamiltonian Q. 

At this point we would like to generate and mea- 
sure the relative phase between the ground and excited 
states. To create the first spin chain excitation, |ejt ), 
in the optical lattice setup we employ two-photon Bragg 
spectroscopy 0. The use of a two-photon process al- 
lows for an accurate energy resolution of the state we 
want to excite, which is crucial considering that the 
energy gap, Afc , is small compared to the tunneling 
and collisional couplings [l9|. For example, an ampli- 
tude modulation of the axial lattice potential, V x (x, t) = 
[Vx,o + A mo< i sin(27rz/ mo di)] sin 2 (fca;), can be employed |20| 
with a frequency that satisfies the resonant Bragg con- 
dition, i/ mo( i = Afe . This will activate a coherent en- 
ergy population transfer from the ground state to the 
excited one with a transition amplitude proportional to 
the modulation amplitude, A mo ^. Such an experimental 
technique has already been employed to study the atom- 
number excitations in an optical lattice near its critical- 
ity between the superfluid and the Mott insulator regime 
[2£j. To detect the Berry phase of the excited state, one 
can initially create the superposition \g) + |efc ) with an 
appropriately timed lattice modulation with frequency 
y mo d = Afc . The Berry phase procedure is then applied 
to generate a phase difference of n between the ground 
and the excited states giving \g) — |efe ). Applying the 
same modulation procedure as before will finally result 
in the state |e/- ). To distinguish between the states \g) 
and |efc ) one can measure the corresponding spin-spin 
correlators, e.g . by atomic scattering of fast atoms off 
the lattice fl7| . 

To conclude, we have presented a method that theoret- 
ically, as well as experimentally, enables the detection of 
regions of criticality through the geometric phase, with- 
out the need for the system to experience phase transi- 
tions. The latter is experimentally hard to realize as the 
adiabaticity condition breaks down and the state of the 
system is no longer faithfully represented by the ground 
state. The origin of the geometric phase can be ascribed 
to the existence of degeneracy points in the parameter 
space of the Hamiltonian. Hence, a criticality point can 
be detected by performing a looping trajectory around it 
and detecting whether or not a non-zero Berry phase has 
been generated. For the case of the XY model the topo- 
logical nature of the resulting phase pinned to the value, 
ip eg « 7T, is revealed by its resilience with respect to small 
deformations of the loop. Topological phases arc inher- 
ently resilient against control errors, a property that can 
be proved to be of a great advantage when considering 
many-body systems. Such a study can be theoretically 



performed on any system which can be analytically elab- 
orated such as the case of the cluster Hamiltonian 
or exploited numerically when analytic solutions are not 
known. The generalization of these results to a wide vari- 
ety of critical phenomena and their relation to the critical 
exponents is a promising and challenging question which 
deserves extensive future investigation. 

Acknowledgements. We would like to thank E. Rico 
for inspiring discussions. This work was supported by 
the EU TopQIP network and the Royal Society. 



[1] M. V. Berry, Proc. Roy. Soc. A 329, 45 (1984). 

[2] A. Shapere and F. Wilczek, "Geometric phases in 

physics.", World Scientific, Singapore, (1989). 
[3] A. Bohm, A. Mostafazadeh, H. Koizumi, Q. Niu, and J. 

Zwanziger, The Geometric Phase in Quantum Systems, 

Ed. Springer (2003). 
[4] V. E. Korepin and A. C. T. Wu, Int. J. Mod. Phys. B 5, 

497 (1991). 

[5] Q. Niu, X. Wang, L. Kleinman, W.-M. Liu, D. M. C. 
Nicholson, and G. M. Stocks, Phys. Rev. Lett. 83, 207 
(1999). 

[6] G. Schiitz, Phys. Rev. E 49, 2461 (1994). 

[7] P. Bruno, "Non-quantized Dirac monopoles and strings 

in the Berry phase of anisotropic spin systems", 

cond-mat/0404616 
[8] J. Wiemer and F. Zhou, "Berry's Phases of Ground 

States of Interacting Spin-One Bosons: Chains of 

Monopoles and Monosegments" , cond-mat/0309312 
[9] D. J. Thouless, Topological Quantum Numbers in Non- 

relativistic Physics, World Scientific, Singapore, London 

(1998). 

[10] S. Sachdev, Quantum Phase Transitions, Cambridge Uni- 
versity Press (1999). 

[11] J. K. Pachos and M. B. Plenio, Phys. Rev. Lett. 93, 
056402 (2004). 

[12] J. K. Pachos and E. Rico, Phys. Rev. A 70, 053620 
(2004). 

[13] A. Kastberg et al, Phys. Rev. Lett. 74, 1542 (1995); G. 

Raithel et al., Phys. Rev. Lett. 81, 3615 (1998). 
[14] A. B. Kuklov and B. V. Svistunov, Phys. Rev. Lett. 90, 

100401 (2003). 

[15] L.-M. Duan, E. Demler and M. D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2003). 
[16] R. Roth and K. Burnett, Phys. Rev. A 67, 031602 (2003). 
[17] A. B. Kuklov, and B. V. Svistunov, Phys. Rev. A 60, 

R769 (1999). 

[18] J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 
Kurn, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 
82, 4569 (1999). 

[19] The energy gap, Ak , is of the order of A, 7. In the optical 
lattice implementation A is manipulated externaly by a 
heavily detuned laser fields, while 7 is of the order of 

Jc / (JaJb)- 

[20] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. 
Esslinger, Phys. Rev. Lett. 92, 130403-1 (2004). 



